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FOREWORD 

This  report  was  prepared  by  Robert  M.  Engle,  Jr.  of  the  Solid  Mechanics 
Branch,  Structures  Division,  Air  Force  Flight  Dynamics  Laboratory.  The  work 
was  conducted  in-house  imder  Project  1467,  "Structural  Analysis  Methods," 

Task  146704,  "Structural  Fatigue  Analysis, "  with  Mr.  Robert  M.  Bader  as 
Project  Engineer. 

This  report  covers  research  conducted  from  July  1969  through  February 
1970. 
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ABSTRACT 

This  report  presents  a  detailed  description  of  a  computer  program  for 
analyzing  crack  propagation  in  cyclic  loaded  structures.  The  program  has  the 
option  of  using  relationships  derived  by  Forman  or  by  Paris  for  crack  growth. 
Provisions  are  made  for  both  surface  flaws  and  '^through  cracks"  as  well  as 
the  transition  from  the  former  to  the  latter.  The  program  utilizes  a  block 
loading  concept  wherein  the  load  is  apjdied  for  a  given  number  of  cycles  rather 
than  applied  from  one  cycle  number  to  another  cycle  number.  Additional 
features  of  the  program  are:  variable  print  interval,  variable  integration 
interval,  and  optional  formats  for  loads  input.  Detailed  input  instructions  and 
an  illustrative  problem  are  presented. 
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SYMBOLS 


Mathematical 

FORTRAN 

Physical 

Symbol 

S3rmbol 

Definition 

0 

A 

Half  crack  length 

»0 

AZERO 

Initial  half  crack  length 

b 

B(2) 

Plate  half -width 

c 

C 

Material  constant 

c 

SMALLC 

Surface  flaw  half-width 

da/dN 

DADN 

Crack  propagation  rate 

KSUBC 

Critical  stress  intensity  factor 

B(3) 

Characteristic  length 

N 

N 

Cycle 

n 

SMALLN 

Material  constant 

R 

R 

Stress  ratio  (  <r^/ 

f 

THICK 

Material  thickness 

BETA(I) 

Individual  correction  factor 

0T 

BETAT 

Total  combined  correction  factor 

AK 

DELTAK 

Stress  intensity  factor  range 

A<r 

TSUBA 

Applied  tensile  stress  range  ( 

^  mb) 

<r 

SIGMA 

Applied  tensile  stress 

Cy, 

SIGMAY 

Material  yield  stress 
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SECTION  I 
INTRODUCTION 

The  total  service  life  of  a  structure  is  often  dependent  upon  the  total 
amount  of  crack  growth  which  can  be  tolerated  prior  to  the  formation  of  the 
critical  size  flaw  or  crack.  An  analysis  which  can  predict  this  growth,  under 
variable  amplitude  loading,  leading  to  the  critical  crack  length  is  a  valuable 
aid  in  establishing  safe  operating  periods  and  inspection  intervals. 

An  automated  procedure  is  presented  in  this  report  which  will  permit  the 
user  to  examine  the  crack  propagation  of  various  flaw  shapes  including  surface 
flaws.  Provision  is  also  made  for  transition  from  a  surface  flaw  to  a  "through 
crack. "  The  computer  program,  CRACKS,  is  written  entirely  in  FORTRAN  IV 
for  the  IBM  7044/7094  Direct  Coupled  System  (DCS).  A  source  listing  is  given 
in  Appendix  I. 
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SECTION  II 

MATHEMATICAL  FORMULATION 

1.  CRACK  PROPAGATION  RATE 

In  the  early  1960’s,  P.  C.  Paris  (Reference  1)  determined  that  the  rate  of 
crack  propagation  under  cyclic  loading  is  primarily  related  to  the  stress- 
intensity -factor  range,  AK.  Paris  proposed  an  exponential  relationship  of  the 
following  form: 

(0 

In  1967,  Forman,  Kearney,  and  Engle  published  a  paper  (Reference  2‘ ,  in 
which  Paris’  equation  was  modified  to  take  into  account  the  effects  of  load  ratio, 
R,  and  crack  growth  instability  as  AK  approaches  K  .  These  modifications  led 
to  a  relationship  of  the  following  form: 

^  CcIAk)"'' 

_1S_  - L -  (2) 

dN  (l-R)  Kc  -Ak 

Both  of  these  relationships  have  proved  useful  in  crack-propagation 
analysis  and  hence  provision  is  made  in  CRACKS  for  both. 

2.  STRESS  INTENSITY  FACTOR 

The  basic  unit  of  fracture  mechanics  is  the  stress  intensity  factor,  K.  For 
crack  propagation  analysis,  the  applied  crack  tip  stress  intensity  factor,  K, 
must  be  less  than  the  material’s  toughness  (K  )  or  fracture  occurs.  This  ap- 
plied  crack  tip  stress  intensity  factor,  K,  is  a  function  of  geometry  and  type  of 
loading.  For  a  central  crack  in  an  infinite  width  plate,  the  stress  intensity 
factor  may  be  written  as  follows: 

K  I  <r  V^TT  0  (3) 


This  equation  will  take  different  forms  based  upon  the  geometry  and  the 
loading.  For  many  cases,  however,  these  effects  may  be  treated  as  modifiers 
or  correction  factors  to  Equation  3.  Thus,  a  more  general  form  would  be: 
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These  correction  factors  will  be  described  in  more  detail  in  the  following 
section. 

Some  investigators  (Reference  3)  have  modified  Equation  4  by  removing  the 
factor  TT  from  under  the  radical  giving: 

K  =  cr  a  /3-j-  (5) 

Equation  4  and  Equation  5  are  both  prevalent  in  the  literature  and  are  in¬ 
cluded  in  the  computer  program. 

The  stress-intensity -factor  range,  AK,  is  defined  as: 

Ak  =  K^ax-^mir. 

Substituting  Equation  4  into  this  relation  gives: 

Ak  s  A  O'  -s/fo  /3t 

Similarly,  substituting  Equation  5  will  yield: 

Ak  =  Ac  /S-j-  (7) 

3.  CORRECTION  FACTORS 

Equations  6  and  7  represent  stress-intensity-factor  ranges  for  a  centrally 
cracked  infinite  panel  if  is  unity.  For  other  geometries,  /3,j,  must  be 
modified.  For  various  combinations  of  geometries,  will  become  combina¬ 
tions  of  different  /8 .  which  will  account  for  these  separate  effects.  For  example, 
in  the  program,  /Sg  corrects  for  finite  width  and  can  correct  for  a  crack 
emanating  from  a  circular  hole.  Hence,  for  a  crack  emanating  from  a  hole  in  a 
finite  width  panel,  )8  ^  would  be  the  product  of  /S  ^  and  13^.  The  program  pro¬ 
vides  for  up  to  ten  P .  of  which  only  four  are  active  at  the  present  time.  Thus, 
in  general, 

Pr  --  TT/3,  (i  =  I,  10)  (8) 
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The  four  active  correction  factors  in  the  program  at  the  present  time 
are  explained  below: 

-  CONSTANT  MULTIPLIER 

This  provides  the  analyst  with  the  capability  to  scale  loads 
or  modify  A  k  by  a  constant  factor 

y3-  -  FINITE  WIDTH  TANGENT  FUNCTION 

This  corrects  for  a  finite  width  plate.  The  form  of  this 
correction  is  (Reference  4) 


where  "a"  and  ••b'*  are  as  shown  in  Figure  1. 


<r 


Figure  1.  Griffith  Crack  in  a  Finite  Width  Plate 

-  -  TABULAR  CORRECTION  FACTOR 

o 

This  permits  the  analyst  to  apply  correction  factors  which 
appear  in  the  literature  as  discrete  data.  The  form  of  this 
correction  is 


=  f  io/J  ) 

O  C 


(10) 
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An  example  is  the  crack  emanating  from  a  circular  hole  (Ref¬ 
erence  4).  In  this  case,  "a'*  and  '*  '*  are  as  shown  in  Figure 

0 

2. 


cr 


/S 


4 


where 


ELLIPTICAL  SURFACE  FLAW  CORRECTION 


The  e^^ression  for  the  stress  intensify  factor  developed  by 
Irwin  (Reference  S)  is  given  as 


1. 1  <r  V IT 0 


K  s 


[<I>^ -0.212  (<r/<7-y,)^] 


2  11/2 


(II) 


The  geometry  of  the  surface  flaw  is  defined  in  Figure  3. 


Figure  3.  Surface  Flaw  Geometry 


As  the  surface  flaw  approaches  the  back  surface  of  the  material,  a  magnifica¬ 
tion  of  the  stress  intensity  factor  takes  place.  This  is  accounted  for  by  a  magni¬ 
fication  factor,  M|^,  which  is  a  function  of  both  a/t  and  a/2c.  The  magnification 
factor  used  in  CRACKS  has  been  obtained  from  Reference  6  and  is  included  in 
Figure  4  for  various  flaw  shapes  (a/2c  values).  Although  this  data  has  been 
derived  for  aluminum,  it  is  in  close  agreement  with  the  results  of  Kobayashi 
(Reference  7)  and  Smith  (Reference  8)  for  general  applications.  Hence,  Equa¬ 
tion  11  may  be  written  as 


where 


K  =  l.l  M|j  <r  ira  /  Q 


(12 


-2  f  <r  \2 

Q  t  (I>  -0.  212  I - ] 


(13 
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Figure  4,  Magnification  Factor  Curves 
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Since  Q  is  a  function  of  er ,  it  is  not  convenient  to  develop  an  equation  of  the 

form  of  Equation  6.  It  is  instead  more  convenient  to  obtain  K  and  K  . 

^  max  min 

and  thus  obtain  AK  .  So,  from  Equations  12  and  13  we  obtain 


•'max  *  ••  -J  ir  a  I 


(14  a) 
(14  b) 


From  these  equations  we  see  that  the  obvious  expression  for  $^i8 

04  <  i-i  Mk  (15) 

and  we  can  then  write 


Ak 


'min 


in 


(16) 


The  translation  from  a  surface  flaw  to  a  through  crack  is  chosen  to  be  the 
point  when  the  plastic  zone  reaches  the  back  face  of  the  material.  The  value  of 
"a"  for  which  this  occurs  is  given  as 


where  is  defined  by  Equation  14a.  At  this  point  an  effective  through-crack 
length  is  calculated  and  the  program  then  continues,  now  using  Equation  6  for 
Ak,  with  0^  set  to  imlty. 
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SECTION  III 

COMPUTER  PROGRAM  FOR  CRACK  PROPAGATION 
ANALYSIS  (CRACKS) 

The  program  described  below  was  written  in  FORTRAN  IV  for  the  IBM 
7044-709411  Direct  Coupled  System.  The  program  consists  of  seven  subpro¬ 
grams,  each  of  which  has  a  specific  task  to  perform.  These  subprograms  and 
their  functions  are: 

CRACKS  -  reads  in  data,  sets  up  calculations,  and  prints  the  results. 

F  -  evaluates  crack  propagation  rate,  da/dN. 

RKIDES  -  variable-step  Runge-Kutta  integration  routine  which  integrates 
da/dN  over  each  load  block. 

TBLKUP  -  linear  interpolation  scheme  for  use  with 

ELIP2  -  routine  to  evaluate  the  complete  elliptic  integral  of  the  second 
kind  to  calculate  <|)  for  use  in  Equation  12. 

TIFANY  -  block  data  subroutine  containing  data  for  Tiffany's  Mj^  curves 
as  a  function  of  a/t  and  a2c. 

TRP2  -  parabolic  interpolation  routine  for  a  function  of  two  variables  which 
is  used  to  determine  Mj^  for  Equation  13. 

A  simplified  flow  chart  depicting  the  transfer  of  information  from  the  sub¬ 
programs  discussed  above  is  given  in  Figure  5. 

I.  DESCRIPTION 

The  CRACKS  computer  program  integrates  the  crack-propagation-rate 
equation  to  obtain  crack  growth  versus  cycles.  The  program  provides  options 
for  the  two  prevalent  forms  of  this  equation  (Equations  1  and  2).  Many  crack 
geometries  may  be  modeled  using  the  correction  factors  discussed  in  Section 

II,  paragraph  3.  A  transition  from  a  surface  flaw  to  a  through  crack  is  pro¬ 
vided  (Equation  17).  As  a  convenience,  the  program  makes  provisions  for  two 
forms  of  spectrum  input,  maximum  and  minimum  stresses,  or  stress  range 
and  load  ratio.  The  program  also  has  the  capability  to  run  multiple  problems 
merely  by  loading  data  decks  in  sequence. 
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Figure  5.  Information  Transfer  Flow  Chart 
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2.  INPUT  INSTRUCTIONS 

The  input  to  the  program  consists  of  four  basic  sections:  analysis  selection, 
material  properties,  geometry  definition,  and  loading.  In  addition,  there  are 
controls  on  print  interval,  repetition  of  load  spectrum,  and  descriptive  data. 
Detailed  instructions  on  inputting  the  data  cards  are  given  below  in  the  follow¬ 
ing  manner: 

a.  Card  number  and  contents 

b.  FORTRAN  name  for  each  variable 

c.  Format  of  card 

d.  Description  of  each  variable 
Card  1  Descriptive  Information 

TITLE 

FORMAT  (14A6) 

TITLE  -  any  information  to  define  the  problem  for  the  user. 

Card  2  Analysis  Selection,  Material  Description 

EQN,  MATID 
FORMAT  (A6,  4X,  7A6) 

EQN  -  alphameric  indicator  governing  choice  of  analysis. 

The  choices  are:  FORMAN  -  use  Equations  2  and  6 
PARIS  -  use  Equations  1  and  6 
NASAF  -  use  Equations  2  and  7 
NASAP  -  use  Equations  1  and  7 

MATID  -  alphanumeric  identification  to  label  material  (may  be  left 
blank) 

Card  3  Material  Properties 

C,  SMALLN,  KSUBC,  SIGMAY 
FORMAT  (4E10.0) 

C  -  material  constant  in  Equations  1  and  2 
SMALLN  -  exponent  in  Equations  1  and  2 
KSUBC  -  fracture  toughness  of  material 
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SIGMAY  -  yield  stress  of  material,  need  only  be  input  when  running 
a  surface  flaw  analysis. 

Card  4  Initial  and  Allowable  Crack  Lengths 

AZERO,  AMAX 
FORMAT  (2E10.0) 

AZERO  -  initial  half  crack  length 

AMAX  -  allowable  half  crack  length.  If  zero,  program  assumes 
infinite  allowable  length. 

Card  5  Correction  Factor  Information 

BETAL,  I,  BI,  BU 
FORMAT  (A4,  II,  5X.  2E10.0) 

BETAL  -  label  for  correction  factors.  Alphameric  characters,  "BETA'* 
I  -  indicator  for  correction  factor  selection.  For  permissible  values 
see  Table  I 

BI  -  variable  for  correction  factor  (see  Table  I) 

BlI  -  secondary  variable  for  correction  factor,  (see  Table  I) 

If  I  equals  three,  a  table  of  points  for  the  correction  factor  follows  the  BETAS 
card  immediately.  The  format  of  this  card  is  as  follows: 

AOVERB,  BETATB 
FORMAT  (2E10.0) 

AOVERB  -  ratio  of  crack  length  to  characteristic  length 
BETATB  -tabular  value  of  /9g  corresponding  to  this  ratio 

Again,  note  that  the  BETAO  card  input  must  always  be  present  and  must  be  the 
last  BETA  card  if  others  are  present. 

Card  6  Initial  Cycle,  Number  of  Applications 

NZERO,  NFLITE 
FORMAT  (ElO.  0,  IIO) 

NZERO  -  number  of  the  initial  cycle  for  this  computer  run. 

NZERO  is  the  cycle  corresponding  to  AZERO. 

NFLITE  -  this  defines  the  number  of  times  the  input  load  spectrum 
will  be  applied. 
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TABLE  I 

CORRECTION  FACTORS 


n 

Corraetlon  Foetor 

BI 

BIX 

Constant  ( ) 

Constont  volua 

— 

2 

Finitt  Width  (  ^2  ) 

Plot*  holf  width 

— 

3 

Tobulor  (/3^  ] 

Chorocttristic  langth 

Numbar  of 
points  in 
tobl* 

4 

Surfac*  flaw  (^84) 

Flow  holf  width 

Matariol 

thicknass 

5-9 

Inoporativt  at  prasant 

— 

— 

0 

End  of  eorraetions 
(Mutt  olways  b«  prastnt ) 

— 

Card  7 


Load  Format  Selector 
LOADS 

FORMAT  (A5) 

LOADS  -  alphameric  indicator  governing  choice  of  input  load  format. 

»'SIGMA"  -  input  or  and  o-  . 

max  min 

"RANGE"  -  input  A<r  and  R 


Card  8  Loads  and  Print  Controls 


LABEL, 


SIGMAX  \  /sIGMInX 

or  ) ,  (  or  )  ,  CYCLES,  NINT,  NPRINT 
DELTAT  /  \  R  / 


FORMAT  -  (A5,  5E10,0) 

LABEL  -  any  five  characters  to  identify  load  block.  May  be  left 
blank.  After  the  last  load  card,  a  card  with  "END"  in 
LABEL  is  required. 

SIGMAX  -  maximum  applied  tensile  stress 
DELTAT  -  applied  tensile  stress  range 
SIGMIN  -  minimum  applied  tensile  stress 
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R  -  stress  ratio 

CYCLES  -  number  of  cycles  in  this  load  block 
NINT  -  integration  interval  for  this  load  block.  If  NINT  is  zero, 
the  program  sets  NINT  equal  to  CYCLES. 

NPRINT  -  print  interval  for  this  load  block.  If  NPRINT  is  zero, 
the  program  sets  NPRINT  equal  to  NINT. 

3.  RESTRICTIONS 

Certain  restrictions  and  limitations  must  be  recognized  or  the  capacity  of 
the  program  will  be  exceeded.  In  general,  violation  of  these  restrictions  will 
not  result  in  termination  of  the  computer  run.  Hence,  stacked  problems  may 
be  salvaged  even  though  an  error  occurs  in  ohe  of  the  first  few  data  decks. 

a.  Program  Capacity 

The  program  will  not  accept  more  than  1000  load  cards  or  100  cards 
for  the  tabular  correction  factor.  This  is  an  arbitrary  choice  of  numbers  and 
may  be  changed  by  modifying  the  appropriate  DIMENSION  statements  within 
the  source  deck. 

b.  Print  Interval 

Because  of  the  variable  step  nature  of  the  Runge-Kutta  integration 
scheme,  it  is  inconvenient  to  generate  printout  at  intervals  less  than  the  inte¬ 
gration  interval.  Hence,  NPRINT  is  always  greater  than  or  equal  to  NINT. 

c.  Negative  Loads 

The  theory  developed  in  References  1  and  2  does  not  permit  negative 
loading.  This  leads  to  the  requirement  that  o-  .  A  cr  »  and  R  must  always 

lU&X 

be  greater  than  or  equal  to  zero. 

4.  OUTPUT 

The  output  generated  by  CRACIffi  is,  for  the  most  part,  self-explanatory. 
The  first  page  of  output  consists  of  a  display  of  the  input  data  which  define  the 
problem.  The  next  section  of  output  consists  of  a  printout  of  the  input  load 
spectrum.  If  the  spectrum  is  input  in  terms  of  maximum  and  minimum  stresses, 
these  are  then  changed  within  the  program  to  stress  ranges  and  stress  ratio  and 
these  are  also  printed  out. 
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After  the  problem  has  been  completely  specified,  the  program  begins 
calculating  the  increase  in  crack  length  by  integrating  either  Equation  1  or 
Equation  2  over  each  load  block  and  printing  results  at  each  print  interval  re¬ 
quested.  This  output  consists  of  values  for  the  most  pertinent  parameters 
which  are  valid  at  the  cycle  printed  out.  If  the  surface  flaw  correction  is  not 
used,  and  Q  from  Equation  14  are  set  equal  to  unity. 

In  the  course  of  the  calculations,  the  denominator  of  Equation  2  is  moni¬ 
tored  to  determine  the  cycle  at  which  it  becomes  negative.  Since  it  is  possible 
for  the  denominator  to  go  negative  during  an  integration  interval,  the  progiram 
sets  an  indicator  in  the  integration  routine  and  retiims  to  the  CRACKS  program. 
Here  the  variables  are  reset  to  the  values  immediately  before  the  denominator 
became  negative.  The  program  then  uses  these  values  as  starting  values  and 
reduces  the  integration  interval.  The  program  then  proceeds  as  before  until 
the  denominator  again  goes  negative.  This  process  is  repeated  until  the  inte¬ 
gration  interval  is  reduced  to  one  cycle.  When  the  denominator  goes  negative 
for  this  integration  interval,  the  cycle  number  at  the  onset  of  instability  is 
available.  The  final  values  are  then  printed  out  and  identified.  A  similar  pro¬ 
cedure  is  followed  when  using  Equation  1  except  that  the  instability  criterion 
becomes  greater  than  nine-tenths  K  . 

IXlaX  C 
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SECTION  IV 
SUMMARY 

An  automated  analysis  for  determining  crack  growth  in  cyclic  loaded  struc¬ 
tures  under  variable  amplitude  loading  has  been  described.  In  Appendix  II  rn 
illustrative  problem  has  been  presented  and  the  results  compared  with  published 
data.  The  results  show  the  program  in  close  agreement  with  the  analytical  re¬ 
sults  and  conservative  in  comparison  with  the  ejqperimental  results. 
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IBFTC  CRACKS  Mq4,XR7,0ECK  CRACOOOO 

CRACOC'Ol 

CRACKS  CRAC00C2 

CRAC0003 

A  FORTRAN  IV  COMPUTER  PROGRAM  FOR  CRACK  PROPAGATION  ANALYSIS.  CRACOOC4 

CRAC0305 

THIS  PROGRAM  IS  RASED  ON  THEORY  DEVELOPED  IN  A  PAPER  APPEARING  IN  CRAC0J06 


THE  'JOURNAL  OH  BASIC  ENGINEERING*  SEPTEMBER  1967, PAGES  459-464, 
ENTITLED 

•NUMERICAL  ANALYSIS  OF  CRACK  PROPAGATION  IN  CYCLIC-LOADED 
STRUCTUR£S',R.G.FGRMAN,V.E.KEARNEY,R.M. ENGLE 

THE  PROGRAM  ACCEPTS  EITHER  FORMAN'S  EQUATION  OR  PARIS'  EQUATION 

FORMAN'S  EQUATION 

DA/ON  »  CIDELTA  K ) ♦*N/ ( ( 1-R IKSUBC-DELTA  K) 

PARIS'  equation 
OA/ON  >  CIOELTA  K}**N 
WHERE 

DELTA  K  «  (DELTA  T)  *  SQRT(PI*A)  •  BETA 

OR  -  FOR  NASA  DATA  - 

DELTA  K  «  (DELTA  Tl  *  SQRT(A)  *  BETA 

AND 

DELTA  T  >(SIGMA  MAX)-(SIGMA  MIN) 

BETA  »  COMBINATION  OF  CORRECTION  FACTORS. 


CRAC00C7 

CRACOOO0 

CRAC0CO9 

CRACOOlO 

CRACOCll 

CRAC0C12 

CRACOOn 

CRAC0Ci4 

CRAC0015 

CRAC0016 

CRAC0017 

CRAC0018 

CRAC00i9 

CRAC0020 

CRAC0021 

CRACOD22 

CRAC0023 

CRAC0024 

CRAC0;'25 

CRAC0026 

CRAC0027 

CRAC0026 

CRAC0C29 

CRAC0030 

CRAC0C31 

CRAC0032 

CRACD033 

CRAC0C34 


INPUT 

CRAC0035 

CRAC0C36 

COLUMNS 

FORMAT 

NAME 

DESCRIPTION 

CRAC0C37 

CRAC0U38 

1-80 

14A6 

TITLE 

ANY  DESCRIPTIVE  INFORMATION 

CKAC0039 

CRAC004C 

1-6 

Ao 

EQN 

•PARIS  '-USE  PARIS'  EQUATION. 

CRAC0vi41 

•FORMAN'-USE  FORMAN'S  EQUATION. 

CRAC0042 

•NASAP'-USE  PARIS'  EQUATION 

CRAC0C43 

WITH  NASA  DELIA  K. 

CRAC0v;44 

•NASAF'-USE  FORMAN'S  EQUATION 

CRAC0:45 

WITH  NASA  DELTA  K. 

CRAC0C46 

7-lD 

4X 

SKIP  4  SPACES 

CRAC0C47 

CRAC0048 

11-52 

7A6 

MATIO 

MATERIAL  ID. ANY  INFORMATION 

CRAC0049 

CRAC005C 

1-lC 

4E10.0 

C 

MATERIAL  CONSTANT  IN  EQUATIONS 

CRAC0C51 

11-20 

smalln 

EXPONENT  IN  EQUATIONS 

CRAC00b2 

21-30 

KSUBC 

CRITICAL  STRESS  INTENSITY  F ACTORCRACOC 53 

31-40 

SIGHAY 

YIELD  STRESS 

CRAC0054 

CRACOObS 

1-10 

2E10.0 

A2EK0 

INITIAL  HALF  CRACK  LENGTH 

CRAC0Cb6 

11-20 

AMAX 

MAXIMUM  CRACK  LENGTH  ALLOWED 

CRACCObT 

CRACOCbH 

1-4 

A4 

8ETAL 

•BETA'-LABEL  FOR  CORRECTION 

CRAC0Cb9 

FACTOR  CARDS 

CRACOOfcO 

5 

11 

I 

CODE  FOR  CORRECTION  FACTORS 

CRAC0061 
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c 

6-lC 

bX 

FIVE  BLANKS 

CRAC0062 

c 

11-20 

2E10.0 

PI 

VARIABLE 

NEEDED 

FOR 

CORRECTION 

CRAC0063 

c 

FACTOR 

CRAC006A 

c 

21-30 

Bit 

VARIABLE 

NEEDED 

FOR 

CORRECTION 

CRAC0C65 

c 

FACTOR 

CKAC0066 

c 

CRAC0067 

C*****  THE  VARIABLES, BI  AND  BII.TAKE  ON  DIFFERENT  MEANINGS  FOR  EACH  BETACRAC0068 
C  CRAC0069 
C  I  CORRECTION  Bl  BlI  CRAC0C70 
C  CRAC0071 
C  1  CONSTANT  CONSTANT  0.  CRAC0072 
C  2  FINITE  WIDTH  PLATE  WIDTH  0.  CRAC0073 

TABULAR 


CONSTANT  0.  CRAC0072 
PLATE  WIDTH  0.  CRAC0073 
CHARACTERISTIC  NUMBER  OF  POINTS  CRAC0074 


LENGTH 


IN  TABLE(NPTS) 


Ct****  THERE  WILL  BE'NPTS’CAROS  TO  DEFINE  THE  TABLE. 
C 

C  1-10  2E10.0  AOVERB  A/8  RATIOS  F 

C  11-20  BETATB  CORRESPONDIN 


A/B  RATIOS  FOR  TABLE 
CORRESPONDING  'BETA' 


VALUES 


SURFACE  FLAW 
END  OF  CORRECTIONS 


FLAW 

HALF-WIDTH 


MATERIAL 

THICKNESS 


C*****  THE  ‘BETAO’  CARO  MUST  ALWAYS  BE  PRESENT. IF  THERE  ARE  SEVERAL 


CRAC0075 

CRAC0C76 

CRAC0077 

CRAC0078 

CRAC0079 

CRACOOBO 

CRAC0081 

CRAC0082 

CRAC0083 

CRAC0084 

CRAC0C85 

CRAC0086 


CORRECTION  FACTORS,IT  MUST  BE  THE  LAST* BETA«CARD. IF  NO  CORRECTIONCRAC0087 


FACTOR  CARDS  ARE  USED, ‘BETAO*  MUST  STILL  BE  PRESENT. 


CRAC0088 

CRAC0089 

CRAC0090 


C  6 

1-10 

ElO.O 

NZERO 

INITIAL  CYCLE  NUMBER 

CRAC0091 

C 

11-20 

lie 

NFLITE 

REPEAT  SPECTRUM  ‘NFLITE*  TIMES. 

CRAC0092 

C 

CRAC0093 

C  7 

1-5 

A5 

LOADS 

LABEL  FOR  FORM  OF  INPUT  SPFCTRUMCRAC0094 

C 

*SIGMA*-READ  CAROIS)  9. 

CRAC0095 

C 

*RANGE'-READ  CARDCS)  10. 

CRAC0096 

C 

CRAC0097 

C  8 

1-5 

A5 

LABEL 

FIVE  CHARACTERS 

CRAC0098 

C 

6-15 

5E10.0 

SIGMAX 

MAXIMUM  STRESS 

CRAC0099 

C 

16-25 

SIGMIN 

MINIMUM  STRESS 

CRACOlOO 

C 

26-35 

CYCLES 

NO.CYCLES  IN  LOAD  BLOCK 

CRACOlOl 

C 

36-A5 

NINT 

INTEGRATION  INTERVAL  IN  CYCLES 

CRAC01C2 

C 

IF  NINT»0,SET  NINT»CYCLES. 

CRAC01C3 

C 

A6-55 

NPRINT 

PRINT  INTERVAL  IN  CYCLES 

CRAC0104 

C 

IF  NPKINT»0,SET  NPRINT«NINT. 

CRAC0105 

A 

MAXIMUM  OF 

1000  LOAD 

CAROS  IS 

PERMITTED. 

CKAC0106 

GO 

TO  CARO  10 

r. 

CRAC0107 

C  9 

1-5 

A5 

LABEL 

FIVE  CHARACTERS 

CRAC0108 

C 

6-15 

5E10.0 

OELTAT 

DELTA  SIGMAISIGMAX-SIGMIN) 

CRAC0109 

C 

16-25 

R 

STRESS  RATIOISIGMIN/SIGMAX) 

CRACOllO 

C 

26-35 

CYCLES 

NO.CYCLES  IN  LOAD  BLOCK 

CRACOlll 

C 

36-45 

NINT 

INTEGRATION  INTERVAL  IN  CYCLES 

CRAC0112 

C 

IF  NINT»0,SET  NINT«CYCLES. 

CRAC0113 

C 

46-55 

NPRINT 

PRINT  INTERVAL  IN  CYCLES 

CRAC0il4 

C 

IF  NPRINT»0,SET  NPRINT=NINT. 

CRAC0I15 

A 

MAXIMUM  OF 

1000  LOAD 

CAROS  IS 

PERMITTED. 

CRAC0116 

c 

CRAC0U7 

C  10 

1-5 

A3 

END 

*END'-TERMINATES  PROBLEM  INPUT. 

CRAC0118 

INTEGER  BLK.ENO 

REAL  NZERO,NINT,NPRINT,NFINAL,NN,NOEL,NSTRr,NWRITE 
REAL  NASAF.NASAP 
REAL  KSUBC,KMAX,MSUBK 


CRAC0119 

CRAC0120 

CRAC0121 

CRAC0122 

CRAC0123 
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COHKON  /DATA/  B(9) ,C,SHALLNtOELTAKtKHAX»KSUBC»>NPTSt JJtSIGHAY»EQN, 
1  IPRC«rHICKtAZEftatAMAXtQ«BETA(<)),lND(9) tPHI,RATIO*SHALLC 

COMHON/CORFAC/TSUBA( 1000) tK(l 000) »AOVERB( IOC) tBETATBC 100) tMSUBK 
COMMON/STOP/ ISrUPt I STOPFf IN, Al, 0X1 tIDX 

DIMENSION  S IGKAXI 1000), SIGMINI 1060), CYCLES (1000), NFINALI 1000), 

1  NINT(1000),NPRINT(100C),TITLE(14) 

DIMENSION  LAB(1000),HAri0(7) 

DATA  END,L0A02  /SHEND  •6HRANGE  / 

DATA  PARIS, FORMAN  /6HPARIS  ,6HFOKMAN/ 

DATA  NASAP,NASAF  /6HNASAP  ,6HNASAF  / 

READ  PROBLEM  SPECIFICATIONS 

10  REA0(5,330)  TITLE 
WRirE(6,400)  TITLE 
KEA0(5,340)  EQN.MATIO 
K£A0(5,3;>0>  C,SMALLN,KSUBC,SIGMAY 
WRITE (6, 410)  MATIO,C,SMALLN,KSUBC,SIGMAY 
REA0(5,350)  AZcRO,AMAX 
IF( AMAX.LE.O.)  AMAX*1.0£37 
WKITE(6,420)  AZERO 
1PRC>0 
00  20  1*1,9 
INb( I)«C 
20  BETA(I)*1.0 

IF ( EON. EQ. PARIS. OK. EQN«EU. FORMAN)  WHITE (6,650) 
IF(EQN.cQ.NASAP.OR.EQN.EQ.NASAF  )  NRITE(6,660) 

READ  CORRECTION  FACTOR  SPECIFICATIONS 

30  REA0(5,370)  BETAL, I ,BI ,611 
IF(l.EQ.O)  GO  TO  100 
GO  TO  (40,50,60,70,80,80,80,80,80),! 

40  BETA(I)*6I 

HRITE(6,450)  BcTA(I) 

INO(I)«I 
GO  TO  30 
50  INO(I)«I 
B(I)«BI 

HRITE(6,610)  B(I) 

GO  TO  30 
60  B(I)«OI 
INO(I)«I 
NPTS»fllI+0.5 
HRITE(6,620)  B(I) 

IF(NPTS.LE.O.OR.NPTS.GT.IOO)  GO  TO  90 
READ! 5, 390)  ( AOVERBI I ) ,BETATB( I ) ,I*1,NPTS) 

MRITE(6,460)(AUVEKB( I),BETATB(1),I*1,NPTS) 

GO  TO  30 
70  B(I)«BI 
THICK*B1I 

RATI0»AZERO/(2.04B(I)) 

SMALLK*(B( I )4*2-AZER0**2)/B(I )*42 

CK*1.-SMALLK**2 

CKSQ0»CK4CK 

CALL  CELI2(PHI0,SMALLK,1.0,CKSQ0,.1ER) 

PHI*PHI0 
IPRC»1 
IND( I)*I 

WRIT£(6,630)  B(I), THICK 
GO  TO  30 
80  CONTINUE 


CRAC0124 

CRAC0125 

CRAC0126 

CRAC0127 

CRAC0128 

CRAC0129 

CRAC0130 

CRACC131 

CRAC0132 

CRAC0133 

CRAC0134 

CRAC0135 

CRAC0136 

CKAC0137 

CRAC0138 

CRAC0139 

CRAC0140 

CRAC0141 

CRACC142 

CRAC0143 

CRAC0144 

CRAC0145 

CRAC0146 

CRAC0147 

CRAC0148 

CRAC0149 

CRAC0150 

CRAC0151 

CRAC0152 

CRAC0153 

CRAC0154 

CRAC0155 

CRAC0156 

CRAC0157 

CRAC015B 

CRAC0159 

CRAC0160 

CRAC0161 

CRAC0162 

CRAC0163 

CRAC0164 

CRAC0165 

CRAC0166 

CRAC0167 

CKACQ168 

CRAC0169 

CRAC0170 

CKAC0171 

CRAC0172 

CRAC0173 

CRAC0174 

CRAC0175 

CKAC0I76 

CRAC0177 

CKAC0I78 

CRAC0179 

CRACOieO 

CRAC0181 

CRAC0182 

CRAC0183 

CRAC0ie4 

CRAC0185 
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GO  TC  30 
90  MftITE(6f640) 

GO  TO  IG 

(<EA0  CONTROL  SPECIFICATIONS 

100  KEA0(5t360)  NZEROtNFLITE 
WRITE(6t430)  NZEKU 
WRITE(6,440)  NFLITE 

READ  LOAD  SPECTRUM  SPECIFICATIONS 

REA0(5.330)  LOADS 
8LK=l 

IF(LOADS.EQ.LOAD2)  GO  TO  140 

LOAD  SPECTRUM  IM  TERMS  OF  MAX  AND  MIN  STRESSES 

READ(5i380)  LAB  (BLK )  .SIUHAX(BLK)  .SIUMINIRLK)  .CYCLES  (C-LK) 
1.  NPRINTIBLK) 

INTEGRATION  INTERVAL  CANNOT  EXCEED  BLOCK  SIZE 

IF (NINT(BLK).GT. CYCLES! BLK) )  NlNT(BLK)«CYCLES(BLK) 
IF(NINT(BLK) .EQ.C.I  NINT(BLK)-CYCL£S(BLK) 
NFINAL(RLK)=NZ=RU+CYCLES(BLKJ 
110  BLK=0LK+1 

READ! 5. 380)  LAB ( BLK ). SIGMAX ( BLK ) .SIGMIN! BLK) .CYCLES ( BLK) 
1.  NPRINTIBLK) 

IF(LAB(BLK).EQ.ENO)  GO  TO  120 

INTEGRATION  INTERVAL  CANNOT  EXCEED  BLOCK  SIZE 


CRAC0186 

CRACOIBT 

CRAC0188 

CRAC0189 

CRAC0190 

CRAC0191 

CRAC0192 

CRAC0193 

CRAC0194 

CRAC0195 

CRAC0196 

CRAC0197 

CRAC0198 

CRAC0199 

CRAC0200 

CRAC0201 

CRAC0202 

CRAC0203 

.NINT(BLK)CRAC0204 

CRAC0205 

CRAC0206 

CRAC0207 

CRAC0208 

CRAC0209 

CRAC0210 

CRAC0211 

CRAC0212 

.NINT(0LK)CRAC0213 

CRAC0214 

CRAC0215 

CRAC0216 

CRAC0217 


IF ( N INT (BLK ) .GT. CYCLES (BLK) )  NINl ( BLK ) "CYCLES (BLK ) 

IF( MINT ( BLK ),EQ.O.)  NINT (BLK) "CYCLES (BLK) 
NFINAL(BLK)=NFINAL(BLK-1)+CYCLES(BLK) 

IF(NPRINT(BLK).EQ.O.)  NPRINT(BLK)"NINT(DLK) 

GO  TO  lie 
120  BLK"BLK-1 
V.RITE(6.470) 

LINE=0 

00  130  I»l,BLK 

WRITE (6.483) LAB (I ).I, CYCLES! I ).NFINAL( I ).SIGMAX(I).SIGMIN(I). 
1  NINT(I) 

THEORY  DOES  NOT  RECOGNIZE  EFFECTS  OF  NEGATIVE  LOADING 

tF(SIGMIN(  I).Lr.O.)SIGMIN( I)"0. 

LINE"LINE+1 
IF(LINE.GT.55)  LINE"0 
IF(LINE.EU.O)  WR(TE(6.470) 

CONVERT  STRESStS  TO  RANGE  AND  RATIO 

TSUBA(I)"SIGMAX(I)-SIGMIN(I) 

R(I)"SIGMIN(I)/SIGMAX(I) 

130  CONTINUE 
GO  TO  170 

LOAD  SPECTRUM  IN  TERMS  OF  STRESS  RANGE  AND  STRESS  RATIO 
140  READ! 5. 380)  LAB(BLK) .TSUBA(BLK) .R(BLK) . CYCLES! BLK ) .NINT (BLK) , 


CRAC021B 

CRAC0219 

CRAC0220 

CRACC221 

CRAC0222 

CRAC0223 

CRAC0224 

CRAC0225 

CRAC0226 

CRAC0227 

CRAC0228 

CRAC0229 

CRAC0230 

CRAC0231 

CRAC0232 

CRAC0233 

CRAC0234 

CRAC0235 

CRAC0236 

CRAC0237 

CRAC0238 

CRAC0239 

CRAC0240 

CRAC0241 

CRAC0242 

CRAC0243 

CRAC0244 

CRAC0245 

CRAC0246 

CRAC0247 
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I  NPklNKBLK) 

INTEGRATION  INTERVAL  CANNOT  EXCEED  BLOCK  SIZE 

IF( NINT (BLK).GT. CYCLES (QLKI)  NINT1BLK)«CVCLES(BLK) 
IFININTIBLKI.EQ.O.)  NINT(BLK)>CYCLES(BLK) 

NF INAL ( BLK ) «NZERO*CYCLES( BLK J 
IFINPRINTIBLKI.EO.O.)  NPKINT(BLK)>NINT( BLK) 

150  BLK>BLK*1 

READ! 5. 360)  LABI BLK ), TSURA ( BLK) »R(8LK) tCYCLES (BLK) , NINT ( BLK) t 
1  NPRINT(BLK) 

IF(LAB(BLK).EQ.ENO)  GO  TO  160 

INTEGRATION  INTERVAL  CANNOT  EXCEED  BLOCK  SIZE 

IF(NINT(BLK).GT. CYCLES  I BLK))  NINT(RLK)«CYCLES(BLK) 
IF(MNT(8LK).EQ.O.)  NINTIBLK) •CYCLESIBLK) 

NFINAL ( BLK ) »NF INAL  I BLK-1 ) ♦CYCLES  I BLK) 

IF(NPRINT(BLK).EQ.O.)  NPRINT(BLK)»NINT(OLK) 

GO  TO  ISO 
160  BLK>BLK*1 
170  URITE(6t500) 

LlNE-0 

DO  160  I«ltBLK 

WRlTE(6fA90)LAB(I)iltCYCLES(I)»NFINAL(I)»TSUBA(I)»R(l)«NINT(l) 

THEORY  DOES  NOT  RECOGNIZE  EFFECTS  OF^NEGATIVE  LOADING 

IFIRin.LT.O.)  R(I)«C. 

L1NE«LINE^1 
(FILINE.GT.SS)  LINE>0 
IF(L1NE.EQ40)  HRITE(6i500) 

160  CONTINUE 
IST0P«0 
IST0PF«0 


CRAC0248 

CRAC0249 

CRAC0230 

CRAC0251 

CRAC0252 

CRAC0253 

CKAC0254 

CRAC02S5 

CRAC0256 

CRAC0257 

CRAC0258 

CRAC0259 

CRAC026C 

CRAC0261 

CRAC0262 

CRAC0263 

CRAC0264 

CRAC0265 

CKAC0266 

CRAC0267 

CRAC0268 

CRAC0269 

CRAC0270 

CRAC0271 

CRAC0272 

CRAC0273 

CRACC274 

CRAC0275 

CRAC02r6 

CRAC0277 

CRAC0278 

CRAC0279 

CRAC0280 

CRAC02ei 

CRAC0282 


IOX>0 

OX1«0* 

A-AZERO 

A1«AZER0 

REPEAT  INPUT  SPECTRUM  NFLITE  TIMES 

00  320  NFLT»l.‘iFLITE 
MRlTE(6i510.  iFLTtA 

IFIEQN.EQ.PARIS.OR.EQN.EG.NASAP)  GO  TO  190 
WRITE(6,520) 

GO  TO  200 
190  WRITE(6t530) 

200  IN-NZERO 
NN»NZERO 
LINES>0 

IFIIPRC.EQ.O)  MKITE(6,560) 

IFI IPRC.NE.O)  HRITE<6t57C) 

DO  LOOP  FOR  LOAD  BLOCKS 

DO  310  JJ«I«BLK 
JJJ>I 

DX-NINTIJJ) 

210  IFIJJ.GT.l)  GO  TO  220 
NSTRT-NZERO 
GO  TO  230 


CRAC0263 

CRAC0264 

CRAC028S 

CRAC02B6 

CRAC0287 

CRAC0288 

CRAC0289 

CRAC0290 

CKAC0291 

CRAC0292 

CRAC0293 

CRAC0294 

CRAC0295 

CaAC0296 

CRAC0297 

CRAC0298 

CRAC0299 

CRAC0300 

CRAC0301 

CRAC0302 

CRAC0303 

CRAC0304 

CRAC0305 

CKAC0306 

CRAC03C7 

CRAC0308 

CRAC0309 
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220  NSTRT=NFINAL(JJ-1) 

230  NWRITE*NSTRT*FLOAT<JJJ)*NPRINT( JJ) 

IF{(NN+NINT(JJ)).LE.NFINAL(JJ))  60  TO  250 
NDEL-NN+NINT ( J J )-NF INAL ( J J  > 

DX*  NINT(JJ)-NDEL 
IF(ISTOPF.EO.O)  GO  TO  250 
240  IFCDXl.GT.l.OC)  IST0PF*0 
DX*DX1 
NN«IN 
A*A1 
IST0P*0 
250  CYC=NN 

INTEGRATE  OVER  ONE  INTERVAL 

CALL  RK1DES(CYC|A,0X) 

IF(ISTOPF.EQ.l)  GO  TO  240 
IF(ISTOP.NE.O)  GO  TO  10 
NN*CYC 

IF(NN.NE.NWRITE)  GO  TO  280 

PRINT  RESULTS  AT  EACH  PRINT  INTERVAL 

CALL  F(CYCtAiOAON) 

IF(IPRC.NE.O)  GO  TO  260 

WKirE(6f540)  NN|A»OAONtOELTAKtKHAX,MSUBKtO 
GO  TO  270 

260  MRITE(6«550)  NN» At SMALLC«OAONfOELTAK,KHAXtHSUBK»Q 
270  JJJ»4JJ+l 

LINES«LINES-i-l 

IF(LINES.EQ*50)LINES*0 

IFCLINES.EQ.O.ANO.EQN.Eti. FORMAN)  WRITE(6t580) 
IFILINEStEQ.O.ANO.EQN.EQ.NASAF  )  NRITE(6t3eo) 
IF(LINES.EQ.O.ANO.EgN.EO*PARIS)  MRITE(6t590) 
IF(LINES.EQ.O*ANO.EQN.EQ«NASAP)  WRITE(6t590) 
IF(LINES.EQ.O.ANO.IPRC.Eg.O)  WRITE f6t 560) 
IF(LINES.EQ.O.ANO.IPRC,NE.O)  WRITE (6t 570) 
IF(NN.EQ.NFINAL(JJ) )  GO  T(.  310 
GO  TO  210 

CHECK  FOR  END  OF  LOAD  BLOCK  AND  PRINT  RESULTS 

280  IF(NN.LT.NFINAL(JJ))  GO  TO  210 
CALL  F(CYCtAtOAON) 

IF( IPRC.NE.O)  GO  TO  290 

WR I TE ( 6 1 540 )  NN i A . OADN t  DEL TAK t KMAX , HSUBK t Q 
GO  TO  300 

290  WRITS(6t550)  NNi At SMALLCt OAONtOELTAKtKMAXtHSUBKtQ 
300  LINES>LINES-M 

IF(L1NES.EQ.50)LINES>0 

IF( LINES. EQ.O. AND. EON. EQ.FOKMAN)  WRITE (6t380) 
IFILINES.EQ.O.AND.EQN.EQ.NASAF  )  WRITE(6t580) 
IFtLINES.EQ.C.AND. EON. EQ. PARIS)  WRITE (6 t 590) 
IFtLINES.EQ.O.ANO.EQN.EQ.NASAP)  WRITE! 6,590) 
IF(LINES.EQ.O.AND.IPRC.EO.O)  WRITE! 6*560) 

IF! LINES.EQ.O. AND. IPRC.NE.O)  WRITE! 6*570) 

310  CONTINUE 

GROWTH*A“AZERO 
WRITE!6t600)  NFLT* AtGROWTH 
320  CONTINUE 

LOOK  FOR  ANOTHER  PROBLEM 


CRAC0310 

CRAC0311 

CRAC0312 

CRAC0313 

CRAC0314 

CRAC0315 

CRAC0316 

CRAC0317 

CRAC0318 

CRAC0319 

CRAC0320 

CRAC0321 

CRAC0322 

CRAC0323 

CRAC0324 

CRAC0325 

CRAC0326 

CRAC0327 

CRAC0328 

CRAC0329 

CRAC0330 

CRAC0331 

CRAC0332 

CRAC0333 

CRAC0334 

CRAC0335 

CRAC0336 

CRAC0337 

CRAC0338 

CRAC0339 

CRAC0340 

CRAC0341 

CRAC0342 

CRACC343 

CRAC0344 

CRAC0345 

CRAC0346 

CRAC0347 

CRAC0348 

CRAC0349 

CRAC0350 

CRAC0351 

CRAC0352 

CRAC0353 

CRAC0354 

CRAC0355 

CRAC0356 

CRAC0357 

CRAC0358 

CRAC0359 

CRAC0360 

CRAC0361 

CRAC0362 

CRAC0363 

CRAC0364 

CRAC0365 

CRAC0366 

CRAC0367 

CRAC0368 

CRAC0369 

CRAC0370 

CRAC0371 
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GC  TL  1C 
FOftK/lT{  14A6) 
FCKf'AT(A6,AX,7A6) 
FGKKflT(6ElC.C) 
FOKf'AT(£i:.CiIlC) 
FOKPAKAA,  Il,5X,2Ci:.:) 
FCRKATC A5f StlC.C) 
FCRKAT(2E1C.CJ 
FGHRAK  IH11CX1AA6) 


CRAC0372 

CRAC0373 

CRAC037A 

CR  ACC  3.7.5 _ 

CRAC0376 

CRAC0377 

CRAC0378 

CRAC0379 

CRAC038C 


FORRAT(lhClCX18HRAT6RIAL  CONSTANTS6X7A6/UX3HC  =E16,8»5X8HSMAtLN  =CRAC0381 


lF5.2t5X7HKSLtJC  »=16.8,5X1AHY1ELD  STRESS  =E16.8) 

I  FOi<KAT(lHClCX27H.NITIAL  HALF  CRACK  LEKGTh  «  E16.8) 

!  FCRRATC  lHClJX22HNmAL  CYCLE  NUKOER  «F1I.3) 

:  FCRRAK lH?i;X21hREP£AT  INPUT  SPECTRUKf I5,7H  TIMES  ) 


CRAC0382 

CRAC0383 

CRAC0384 

CRAC0385 


!  FCRMAT(1HCICX3AFOETA(1)-CONSTANT  CORRECTION  FACTOR  /19X9HDETA(1)  sCRAC0386 
1  E16.8)  CRA.C0387.. 

:  FCRPATI 1F013X3FA/H1£X7HUETA(3)//I8X£15.8,5XE15.8) )  CRAC0388 

:  FORMAT! 1H112XAFLCAD11X9HNUMBER  0F11X9HNUMBER  0F18X5HSIGMA20X5HSIGMCRAC0389 
lAlAXllFINTEGRATIC\/13X5hBLCCKlCXlCHCYCLES  PER1CX8HCYCLE  AT19X3HMAXCRACC390 
222X2FMIN16X8HINTERVAL/2&X1CHL0AD  BLOCK10X12HEND  OF  BL0CK59X8H(CYCLCRAC0391 
3cS)//)  CRAC0392 

:  FCRMAT(2XA5,6X15»10XFil.3,l0XF11.3f 15XF1C.2,13XF10.2,IOXF11,3)  .  C8AC0.393 

3  FORMAT(2XA5|6XIS,10XF11.3,10XF11.3tl5XFlC.2.16XF6.3t  10XF11.3I  CRAC039A 
:  FORMAT! lHll2XAhLCADllX9hNUMBER  0F11X9HNUMBER  CFl8X5hDELTA22XlHRl6XCRAC0395 
lUHNTEGRATlCN/12X5h8LCCKlCXlOHCYCLES  PERIOX8HCYCLE  AT19X5HS IGMA39CRAC0396 
2X8hIATERVAL/28XlCHLCAn  liLCCKlOX12HEND  CF  BLLCK59X8HICYCLES)//)  CRAC0397 
3  FCRMAT!lHl5X35hCRACK  LcNCTh  AT  BEGINNING  OF  FLIGHT  I5»3H  ISF10.5/)CRAC0398 
:  FORMAT! lH?ACX5ChCRACK  PKCPAGATICN  ANALYSIS  USING  FORMAN'S  EQUAT I0NCRAC0399 
1  /AAXAAFCA/DN*C#!0ELTA  K)**N/(!l-R)*KSUBC-CELTA  K)  )  CRAC0400 

J  FORMAT! 1h:41X48FCKACK  PROPAGATION  ANALYSIS  USING  PARIS'  EQUATION  /CRAC0401 
1  55X23FOA/DN=C*!OELTA  KJ^ON  )  CRAC0402 

:  FORMAT! lXFll.3,12XF9,5,13XEl5.8,13XFlC.2»l6XFlC.2t4XF5,3t5XF6.3)  CRA00403 

3  FORMAT! lXF11.3,13XF9.5.2XF9.5t2XE15.8,13XFi:.2,l6XFlC.2tAXF5.3,5X  CRA06404 
1  F6.3)  .  CRA.C04Q5 

:  FORMAT! lMC2XbHCYCLE20XlHA22X5HOA/DN21X7HCELTA  K19X5HK  MAX6X5HMSUBK0RA00406 
1  aXlFC//)  CRAC0407 

3  FORMAT! lH02X5HCYCLE2CXlhAlCXlHCllX5HDA/DN21X7HDELTA  K19X5HK  MAX  CRAG0408 
1  6X5FMSUeK8XlFQ//)  CRAC0439 

;  FORMAT! IHIA’-XSCHCRACK  PRCPAGATICN  ANALYSIS  USING  FORMAN'S  EQUAT I0NCRAG041C 

1  /44X44F.DA/CN»G*!DELTA  K  l-R)  ♦KSUBC-DELTA  .K )  ).  _ _ CRACQ411 

3  FORMAT! IH141X48FCRACK  PROPAGATION  ANALYSIS  USING  PARIS'  EQUATION  /CRAG0412 
I  55X23HCA/CN=C»!DELTA  K)**N  }  CRAG0413 

:  FORMAT! lHCi:x3eFCRACK  LENGTH  AT  END  OF  FLIGHT  NUMBER  15, 3H  ISF1C.50RAC0414 
1/  11X21FTCTAL  CKAOK  GRCNTH  IS  FIC.5) 

:  FORMAT! lhClCX3eFReTA!2»-FIMTE  «IDTH  CCRREOTION  FACTOR  / 19X42HBETACRAC0416 
U2)  =  SQRT!2/!Pl*A/e)*TAN(PI*(A/B)/2))  /l9Xl2HhHERE. ....a  =-.Elj6.a.)..CRAC0.417.. 
3  FORMAT! lHCKX32FeETA! 3)-TABULAR  FUNCTION  OF  A/B  /l9Xl2hWHERE....BCRAC0418 
1  =  rlE.S)  C.RAC.0.419  . 

3  FORMAT!  lhClCX39FHETA!4)-SURFACE  CRACK  CORRECTION  FACTOR  /  CRAC0420 

1  19X51FC£TA!4)=1.1/S0RHPHI**2-C.212<‘!DELTA  T/SIGMA  Y)**2)/  CRAC0421 

219X51FNh£RE....PFI=E!K,PI/2),K«! !B**2-A**2)/B**2)  AND  B  *  E16.8  /CRAC0422 

324X62F.E!K,PI/2)  IS  TFE  COMPLETE  ELLIPTIC  INTEGRAL  CF.  THE  SECOND  KICBACQ423... 


4NC/lSX26hMATERIAL  ThICNNF.SS . T  »  E16.8) 

3  FORMAT! 1HCICX33HINPLT  ERROR. NPTS*0, OR  NPTS.GT.lOO  ' 
3  FORMAT!  1HC1:X36FC0RRECTICN  FACTORS  USED  FOR  DELTA  K 
1  K  =!CELTA  T)*SCRT(PI*A)*BETAIi)*BETAt2)4...BETA(I) 
3  FORMAT! 1HC1CX36FCURRECTILN  FACTORS  USED  FOR  DELTA  K 
1  K  »!CELIA  T)*SCRII  A  )*DETA! 1) ♦BETA! 2 ) ♦. . .BETA ( 1) 
END 


CRAC0424 

CRAC0425 

//11X60F0ELTACRAC0426 

»CBAC0427 

//llX60hDELTACRAC0428 

.)CRAC0429 

CRAC043C 
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(IBFTC  DADN  H94tXR7fDECK 

SUBROUTINE  F(NtA,OAON) 

COHNGN  /DATA/  B(9) ,CtSHALLNtOELTAK»KMAX»KSUBCtNPTSt JJ*SIGMAYtEONt 
1  IPKCfTHlCK.AZER0.ANAXtQ«8ErA(9) »IND(9),PHltRAriO»SHALLC 

COHMON/CORFAC/TSUBAdOOO.RI  1000)  tAOVERBI  100)  *BETATB(100)tKSUBK 
COMHON/STOP/lSTOPt ISTUPF. IN«Al«OXlt lOX 
COMMON  /  MKCURV  /  MK(IOOO) 

REAL  MKtKSUBKtHSUBKl 

REAL  NASAF.NASAP 

REAL  KSUBCfNtKMAXfKHAXl 

REAL  KHIN 

REAL  INtINN 

DATA  PI/3. 1415926/ 

DATA  PARIS. FORMAN  /6HPARIS  .6HF0RMAN/ 

DATA  NASAP.NASAF  /6HNASAP  .6HNASAF  / 

BETAT*1.0 

MSUBK*1.C 

0*0.0 

SMALLC*A/(2.0«RATI0) 

DO  70  1*1,9 
J*IND(I) 

IF(J.EO.O)  GO  TO  70 
GO  TO  (60,10.20.40.50.5C.50.50.50),J 
10  IF( IPRC.NE.C)  GO  TO  70 
SMALLK*A/B{ J) 

IFISMALLK.GT.l.O)  GO  TO  320 
BETA(J)=SttRT(2./(PI*SMALLK)*TANIPI*SMALLK/2.)) 

GO  TO  6C 

20  SMALLK*A/e(J) 

30  BETA(J)sr8LKUP(A0VERB,BETATB,NPTS»100,SMALLK) 

GO  TO  60 

40  ATRANS*THICK-( ( (XMAX/SI6MAY)**2)/(2,0*P1 ) ) 

IFLA.GE.ATRANS)  GO  TO  340 

AOVERT  *  A/THICK 

MSU8K  «  TRP2(MK, AOVERT, RATIO.l) 

BETA!  J)*1.14>HSUBK 
GO  TO  60 
50  CONTINUE 

60  BETAT=BETAT*BETA(J) 

70  CONTINUE 

IF!  IPRC.EQ.O)  GO  TO  80 

Q  *  SHAPE  FACTOR  FOR  SURFACE  CRACKS 

S  IGMAX=TSUBAU  J  )  /  ( 1 . 0-R  ( J  J ) ) 

SIGMIN»SIGMAX-rSUOA(JJ) 

0MAX*PHI**2-0.212*(SIGMAX/SIGMAY)**2 
0*0HAX 

0MIN*PHI**2-0.212*(SIGMIN/SIGMAY)**2 
80  IFIEON.EO.NASAP.OR.EQN.EQ.NASAF)  GO  TO  100 

DELTA  K  USING  FORMAN'S  FORM  AND  CONSTANTS 

IFIIPRC.NE.C)  GO  TO  90 
OELTAK*TSUBA(JJ)*SQRT(PI*A)*BETAT 
KMAX*DELTAK/(1.0-K(JJ)) 

GO  TO  120 

90  KMAX=SI0MAX*SQRT(PI#A/QMAX)*8ETAT 
KMINsSIGMINfSGRTIPI+A/QMINI^BETAT 
DELTAK*KMAX-KMIN 
GO  TO  120 


DADNCCOO 

OADNOOOl 

OADN0002 

0ADN0003 

0A0N00C4 

DADN0005 

OADN0006 

DAON0007 

0ADN0008 

0ADN00G9 

DADNOOIO 

OADNOCll 

0ADN0012 

DA0N0013 

0A0N0014 

0ADN0015 

DA0N0016 

0A0N0017 

0ACNC018 

DADN0019 

DADN0020 

DADn6021 

DADN0022 

DADN0023 

0ADN0024 

DADNOw.25 

DA0N0Q26 

0ADN0027 

DADN0028 

DADN0029 

DA0N003C 

DADN0C31 

0ADN0O32 

DADN0C33 

DA0N0034 

0A0N0U35 

0ADN0036 

0ADN0037 

OAON0O3B 

0ACN0039 

DA0N0040 

0ADN0041 

DA0N0C42 

OA()N0043 

DADN0044 

DA0N0045 

DADN0046 

0ADN0C47 

DADN0048 

DADiN0049 

OADN0050 

DADN0051 

0ADN0O52 

DA0N0053 

DADN0054 

DADN0055 

0ADN0C56 

DADN0057 

DADN00S8 

0ADN0059 

DA0N0060 

0A0NC061 
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C 

C 


DELTA  K  USING  HUDSON'S  FORM  AND  CONSTANTS 

100  IF( IPMC.NE.O)  GO  TO  110 

OELTAK«rSU6A(JJ)*SORT(A)*BETAT 

KHAX«0ELTAK/(1.0-R(Jjn 

GO,.  TO  120 

110  kMAX«SIGMAX*SQRT(  A  /OMAXI^OETAT 
KHIN-SIGH1N«SQRT(  A  /QMINI^BETAT 
OELTAK«KHAX-KMIN 

120  IF(A.GT.AMAX)  go  TO  360 

IFtEGN.EQ. PARIS. OR. EON. EO.NASAP)  GO  TO  130 
OENOM« ( 1 .O-R ( J J ) ) PKSURC-OELTAK 
IFIDENOH.LE.J.)  GO  TO  150 
GO  TO  lAO 

130  DELKC*0.9<»KSUBC-KHAX 

IFIOELKC.LE.O.I  GO  TO  15C 
OENOH-l.O 

lAO  OAON* ( C*DELT AK*«SHALLN ) /DENOH 
IN«N 


OADN0062 

DADN0G63 

OAONOC6A 

OADNOQ65 

0ADN0C66 

0A0N0067 

OA0N0068 

DA0N0069 

0ADN0070 

DAUN0071 

0ADN0072 

0A0N0073 

DA0N007A 

OA0NO075 

DADN0C76 

0ADN0077 

OADN0076 

OA0N0079 

OADN0080 

0ADN0081 

DA0N0082 

0ADN0083 

DA0N0084 

0A0N00B5 

OADN0086 

OADN0087 


GO  TO  290  0ADN0081 

150  ISrOP>l  DA0N0082 

IF(ISTOPF.NE.O)  GO  TO  21C  0ADN0083 

tSTGPF«IST0PF4l  DA0N0084 

IF(IOX.NE.O)  GO  TO  180  0A0N00B5 

INN«N  OADN0086 

IFIEGN.EQ. FORMAN. OR. EON. EO.NASAF)  WR1TE(6*160)  INN  OADN0087 

160  FORMAT! IH036H((1>R)KSUBC-OELTAK)  mENT  NEGATIVE  ATF11.3»7H  CYCLES  /DA0N0088 

llH010X37HnEGIN  SEARCH  FOR  MOKE  ACCURATE  VALUES  /11X15HSTART ING  VADA0N0089 
2LUES//)  OADN0090 

IF(EQN.EQ.PARIS.OR.EQN.EO.NASAP)  NR1TE(6»170) INN  DADN0091 

170  FORMATaH037HK  MAX  BECAME  GREATER  THAN  .9*KSUbC  ATF11.3»7H  CYCLES/OADN0092 

11H01CX37HREGIN  SEARCH  FOR  MORE  ACCURATE  VALUES  /I1X15HSTARTING  VAOA0N0093 
2LUE$//)  DADN0094 

GO  TO  200  OADN0095 

180  WRITE(6»19C)  lOX  OA0N0096 

190  FORMAT!  IHOIOXUHVALUES  FOR  OX  «  14*  7H  CYCLES  //)  DADN0097 

200  HKITE!6*220)  IN.AliOAl,OELl*KHAXltMSUBKl*Ql  DADN009B 

GO  TO  240  OADN0099 

210  MKlT£!6f230)  OADNOlOO 

HRITc!6»220l  IN*A1*0A1,0EL1*KMAX1*HSUBK1*Q1  OAONOlOl 

220  FORMAT!2XFll.'i«13XF9.5«13XE15.Btl3XF10.2*16XF10.2t4XF5.3*5XF6.3)  DA0N0102 

230  FORMAT! IH010X41H**44*AVALUES  AT  UNSET  OF  INSTABILITY*****  ///»  OA0N01O3 

1ST0PF«2  0ADN0104 

RETURN  0A0NC105 

240  OELTAN*N-IN  DADN0106 

IFIDELTAN.LE.l.)  GO  TO  210  0ADN0107 

N0EL«AL0G10!0ELTAN)  0A0N0108 

IFINOEL.EQ.C)  GO  TO  250  0ADN0109 

IFIN0EL.GT.6)  GO  TO  300  OAONOllO 

GO  TO  !250.250*230*260»270»260)*NDEL  OAONOlll 

250  0X1>1.  OAON0112 

IOX>OX140.5  DADN0113 

RETURN  0ADN0114 

260  0X1*10.  DAUN0115 

IOX-0X140.5  OA0N0116 

RETURN  DA0N0117 

270  0X1>10C.  0ADN0118 

10X*DX140.5  OADN0119 

RETURN  •'  0A0N0120 

280  0X1-1000.  DA0N0121 

IDX-0X140.5  DADN0122 

RETURN  OADN0123 
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290  A1«A 

OAl-OADN 

0EL1«0ELTAK 

KMAXl-KNAX 

Ql«(l 

MSUBKl-HSUBK 

KETURN 

300  MRITE(6t310i 


0ADN0i2A 

0ADN0125 

DADN0126 

0ADN0127 

DA0N0128 

0A0N0I29 

0A0N0130 

0A0N0I31 


310  F0RMAT(1H06X101HS0METH1NG  IS  RADICALLY  WRDNG.THE  CRITICAL  VALUE  0FDADN0132 


1  N  CAN  ONLY  BE  DETERMINED  WITHIN  ONE  MILLION  CYCLES  ) 
ISTOPF-2 
RETURN 

320  WRITE(6»330) 


0A0N0133 

OADN0134 

DADN0135 

0A0N0I36 


330  FORMAT (1H075HCRACK  LENGTHtAtiS  GREATER  THAN  PLATE  WIDTH»8.THIS  IS  0ADN0137 
IPHYSICALLY  IMPOSSIBLE.  )  OAONOI38 

GO  TO  390  DADN0139 

340  AEFF«8ETAT«*2*ATRANS  OA0NO14O 

URITE(6|350)  AEFF  OA0N0141 

350  FORMAT! lX130(lH4)/35X41HTRANSiTI0N  TO  A  THROUGH  CRACK  OF  LENGTH  0ADN0142 
1  F9.5t7H  INCHES/1X130(1H*I)  OA0N0143 

IN0t4)>0  0A0N0144 

A>AEFF  DADN0145 

IPRC-0  0ADN0146 

GO  TO  80  0ADN0I47 

360  WRITE(6*370I  DADN0148 

370  FORMAT! 1H1130(LH*)/14X104HTHE  CRACK  GROWTH  EXCEEDS  THE  MAXIMUM  ALLOA0N0149 
LOWED  LENGTH.  EXAMINE  THE  PROBLEM  FOR  POSSIBLE  REFORMULATION.  /DADN0150 
21X130!  lH*n  DADN0I51 

WRITE!6t380i  0ADN0I52 

380  FORMAT! IH03X5HCVCLE20XIHA22X5HDA/0N21X7HDELTA  K19X5HK  HAX6X5HMSUBKDADN0153 


I  8X1HQ//)  0A0N0154 

WRtTE!6f220l  IN«AItOAl(OELltKMAXl»MSUBKltQl  0ADN0155 

390  ISTOPsl  0ADN0156 

ISTOPF-2  0ADN0157 

RETURN  DAONOlSe 

END  0ADN0159 
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SIBFTC  RKIOES  M94iXI<7,DECK 

SUBKOUTINE  i<KIOES(X,YtOX) 

COHHON/STOP/ISrOPt ISTOPFt INtAltDXlt lOX 
10  XO«X 
X-X*DX 
H»OX 

20  IF(ABS(H).GT.A8S(X>X0) )  H>X-XO 
30  YO»Y 
HT»H 
XT«XO 

RHAXP-1.E37 
40  YT-VO 

ASSIGN  50  TO  K 
GO  TO  100 
SO  CONTINUE 
60  YP«Y 
70  HT*0.5*H 

ASSIGN  80  TO  K 
GO  TO  100 
80  CONTINUE 
90  YT«Y 
XT-XO^HT 
ASSIGN  150  TO  K 
100  CALI.  F(XT,YT,PO) 

IFi ISTOP.NE.O)  RETURN 
110  Y*YT*0.5*HT*P0 

CALL  F(XT+0.5*HT,YfPll 
IFI ISTOP.NE.O)  RETURN 
120  Y*YT+0.5*HT*P1 

CALL  F(XT*0.5*HT,YfP2) 

IFdSTUP.Nc.G)  RETURN 
130  Y*YT+HT4P2 

CALL  F(XT*HT»YtP3) 

IFI  ISTOP.NE.O)  RETURN 
140  Y*YT+HT*(P0*2.*(Pl+P2)+P3)/6. 

GO  TO  Ki(50»80il30) 

150  KMAX«0. 

160  RMAX»ARAXHRMAX.0.07*ABS((Y-YP)/Y) ) 

IFI (RMAX.GT.l.e-C6).AN0.(RMAX.LT.RMAXP) )  GO  TO  170 
X0»X0»H 

IFIXO.EG.X)  RETURN 

IFI IKMAX.lt. l.E-07).  OR. (RMAX.GT.RMAXP) )  H-H*H 
GO  TO  20 
170  H»HT 
XT»XO 
180  YP«YT 
190  YT»YO 

RHAXP>RMAX 
GO  TO  70 
END 


RKIOOOOO 

RKlOOOOl 

RK100002 

RK100003 

RK100004 

RK100005 

RK100006 

RK100007 

RK1000C8 

RK100009 

RKIDOOIO 

RKlOOOll 

RK100012 

RK100C13 

RK100014 

RK100015 

RK1D0016 

RK100017 

RK100018 

RK1D0019 

KK100020 

RK100021 

RK100C22 

RK100023 

RK100024 

RK1D0025 

RK100026 

KK100027 

RK1D0028 

RK1D0029 

RK1D0030 

RK100031 

KK1D0032 

RK100033 

KK100034 

RK100035 

RK1D0036 

RK100037 

RK100038 

RK 100039 

RK100040 

RK100041 

RK100042 

RK100043 

RK1D0044 

RK100045 

RK100046 

RK100047 

RK1D0048 

RK1D0049 


30 


AFFDL-TR-70-107 


SIBFTC  TBLKUP  M94fXR7iOECK  TBLKOOOO 

FUNCTION  TDLKUP(XtYtNtNHAXtARG)  T8LK0031 

DIMENSION  X(NHAX)fY(NHAX)  TBU0002 

DO  10  I<lfN  TBLK0003 

IF(X(I)-ARG)  10t20f20  T6LK0004 

10  CONTINUE  TBLK0005 

I«N  TBLK0006 

20  IF( 1-1)30,30*40  TBLK0007 

30  1«2  TBLKOOOB 

40  SLOPE»(Y(n-Y(l-l))/(X(I)-X(I-l))  TBLX0009 

TBLKUPsSLOPE*(  ARG-X(  I-l)  )+Y(  I-l)  TBUOOlO 

RETURN  TBLKOOIl 

END  T6LK0012 
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SIBFTC  EL1P2  M94tXK7tOECK 


SUBROUTINE  CEL  12 
PURPOSE 

COMPUTES  THE  GENERALIZED  COMPLETE  ELLIPTIC  INTEGRAL  OF 
SECOND  KINO. 

USAGE 

CALL  CELI2(RES«Ar.«A«B»lER) 

DESCRIPTION  OF  PARAMETERS 
RES  >  RESULT  VALUE 

AK  -  MODULUS  (INPUT I 

A  *  CONSTANT  TERM  IN  NUMERATOR 

8  -  FACTOR  OF  CUAORATIC  TERM  IN  NUMERATOR 

lER  -  RESULTANT  ERROR  CODE  MHERE 

lcR«0  NO  ERROR 

1ER«1  AK  NOT  IN  RANGE  >1  TO  ♦! 

REMARKS 

FOR  AK  •  ♦It-l  THE  RESULT  VALUE  IS  SET  TO  1.E75  IF  B  IS 
POSITIVE*  TO  -1.E75  IF  R  IS  NEGATIVE. 

SPECIAL  CASES  ARE 

K(K)  OBTAINED  WITH  A  •  1*  U  >  1 

E(K)  OBTAINED  WITH  A  •  1*  U  •  CK4>CK  WHERE  CK  IS 

COMPLEMENTARY-  MODULUS. 

6(K}  OBTAINED  WITH  A  >  1*  B  •  0 

0(K)  OBTAINED  WITH  A  «  0*  B  >  1 


WHERE  K*  E«  B*  0  DEFINE  SPECIAL  CASES  OF  THE  GENERALIZED 
COMPLETE  ELLIPTIC  INTEGRAL  OF  SECOND  KINO  IN  THE  USUAL 
NOTATION.  AND  THE  ARGUMENT  K  OF  THESE  FUNCTIONS  MEANS 
THE  MODULUS. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

METHOD 

DEFINITION 

RES» INTEGRAL ( (A*B*T*TI/(SQRT((l+T*T>*(U(CK*T)**2»)*(l*T*T» 
SUMMED  OVER  T  FROM  0  TO  INFINITY). 

EVALUATION 

LANOENS  TRANSFORMATION  IS  USED  FOR  CALCULATION. 

REFERENCE 

R.OULIKSCH,  'NUMERICAL  CALCULATION  OF  ELLIPTIC  INTEGRALS 
AND  ELLIPTIC  FUNCTIONS**  HANDBOOK  SERIES  SPECIAL  FUNCTIONS* 
NUMERISCHE  MATHEMATIK  VOL*  7*  1965*  PP.  78-90. 


SUBROUTINE  CELI2(RES*AK*A*8*IER) 
IER»0 

TEST  RANGE 
CK«AK*AK 

IF(CK-1.)  20*20.10 
10  IER«1 
RETURN 


ELIPOOOO 

ELlPOOOl 

.ELIP00C2 

EL1P0C03 

EL1P0004 

EL1P0005 

ELIP0006 

ELIP0007 

ELIP0008 

EL1P0009 

ELlPOOlO 

ELIPOOll 

ELIP0012 

EL1P0013 

ELIP0014 

EL1P0C15 

EL1P0016 

ELIP0317 

ELIP0018 

ELIP0019 

ELIP0020 

EL1P0021 

EL1P0022 

EL1P0023 

EL1P002A 

ELIPOOZS 

ELIP0026 

EL1P0027 

ELIP0028 

ELIP0029 

ELIP0030 

ELtP0C31 

ELIP0032 

ELIP0033 

EL1P0034 

ELIP0035 

ELIP0036 

ELIP0037 

ELIP0038 

ELIP0039 

EL1P0040 

)ELIP0041 

EL1P0042 

ELIP0043 

ELIP004A 

EL1P0C45 

EL1P0046 

ELIP0047 

ELIP0048 

EL1P0049 

.ELIP0050 

ELIP0051 

ELIP.0052 

EL1P0053 

ELIP0054 

ELIP00S5 

ELIP0056 

ELIP0057 

EL1P0058 

ELIP0059 

ELIPOOOO 

ELIP0C61 
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COMPUTE  COMPLEMENTARY  MODULUS 

20  GE0*SQRT(1.C-CK) 

IF(GEO)  70t30f70 

SET  RESULT  VALUE  >  OVERFLOW 

30  1F(D)  40t60tS0 
AO  RES*-l.E38 
RETURN 

50  RES>1.E38 
RETURN 
60  KES-A 
RETURN 

COMPUTE  INTEGRAL 

70  ARI^l. 

AA>A 

AN«A4B 

w«a 

80  W-W4AAKGE0 
H*W4W 
AA«AN 
AARI*ARI 
AR1>GE04ARI 
AN«W/ARI4AN 

TEST  OF  ACCURACY 

IF(AARI-GE0>l.E-4«AAKn  lOOflOOtOO 
90  GE0«SQRT(GE0*AARn 
GE0>0E04GE0 
GO  TO  80 

100  RES».78539816*AN/ARI 
RETURN 
END 


ELIP0G62 

EL1P0063 

ELIP0064 

ELIP0065 

ELIP0U66 

ELIP0067 

ELIP0068 

ELIP0069 

ELIP0070 

ELIP0071 

ELIP0072 

ELIP0073 

EL1P0074 

ELIP0075 

ELIP0076 

ELIP0077 

EL1P0078 

EL1P0079 

ELIP0080 

ELIP0081 

ELIP0082 

ELIP0083 

EL1P0084 

ELiPOOSS 

ELIP0C86 

ELIP0087 

ELIP0CB8 

EL1P0089 

ELIP0090 

ELIP0091 

ELIP0092 

ELIP0093 

ELIP0094 

ELIP0095 

ELIP0096 

ELIP0097 

EL1P0098 
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*IBFTC  TRf»2  M94,XR7»0ECK 

REAL  FUNCTION  TRP2(TtX»Y»M) 

OIMcNSlGN  T(  1003ltZ(4)«U(6t 

Ll«C 

Xl«X 

Yl«Y 

I»T<l)/lCCC.*l, 

J*AHCO(T(l)»10OC.I^l. 

L«J«M 

Ii»J*3*l 

I2«I*J 

Ml«M 

00  IC  K>lltl2iL 
IF(Xl>r(K))  20i20*lC 
13  CONTINUE 
K-I2+1-J 

20  00  30  L«4«JtMl 

IF  (Y1>T(L))  A0t40t30 
30  CONTINUE 
L*J 

AO  L1«L1*1 

DO  50  MN«li3 
N«L+MN“3 

Ni«K4U*(Ll-'3n*N-l 

0(MNI«T(N) 

50  0(KN«3)*T(N1) 

60  2(Ln«0(A)«(Yl>0(l))«((0(5)*0(4))AIDI2)>D(l))^( 
iYi-o(2n/<o(3)-o(i))«no(6)«o(5n/(o(3)*o(2n 
2>(o(S)-o(An/(o(2)«o(un) 

IF  (Ll-3H0»70»90 
70  00  80  MN«1«3 
0(MN*3)«Z(NN) 

Nl«K4(-JA(MN-3)  J 
80  0(HN)>T(NU 
Ll«4 
Y1«X 

CO  TO  60 
90  TRP2»Z<A» 

RETURN 

ENO 


TRP20000 

TRP20001 

TRP20002 

TRP20003 

TRP2000A 

TRP20005 

TRP20006 

fRP20007 

TRP200C8 

TRP20009 

TRP20010 

tRP20011 

TRP20012 

TRP20013 

TRP20014 

TRP20015 

TRP20016 

TRP20017 

TRP20018 

TRP20019 

TKP20020 

TRP20021 

TRP20022 

TRP20023 

TRP20024 

TRP20025 

TRP20026 

TRP20027 

TRP20028 

TRP20029 

TRP20030 

7RP20031 

TRP20032 

TRP20033 

TRP20034 

TRP20035 

TRP20036 

TRP20037 

TRP20038 

TRP20039 
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$I8FTC  TIfANY  H94tXk7tnECK 
BLOCK  DATA 

COMMON  /MKCURV/  MK(IOOO) 

REAL  MK 

OATA(MK(  n,I>l«fiA)/ 

llC06.«0.05t0.10i0.20tO*30»C.AOt0.50t 

0*0il«00fl*0Ufl*  00  *  1 • 00 1 1*00  y 1 • 00 1 
O.lyl.Oltl.Oltl.OItl.Ol.l.CItl.OOf 
0.2il.03,1.03il.02»1.02tl.01tl.00y 
0.3,1.06.1.06«l.C4»1.03.1.02»UOOt 

C*Af I«12tl*12il*08f 1*05  y 1*02  y 1*  00  y 
0.5.1.22yl.l8yl.l4yl.08yl.03yl.00y 
0«6yLtt34yly30yI«22yl«13yl*06yl*01y 
0«7yl«48yly42fl«31yl-*20yl*08yl*02y 
G.8yl.64yl.57yl.41yl.26yl.l3yl.04y 
0.9yl.77yl.66yl.^0yl.32yl.l8yl.0ay 
l.Of l*84yl.75f 1.59yl.38yl.22yl.l0/ 

END 


TIFAOOCO 

TIFAOOOl 

TIFA00G2 

TIFA0003 

TIFA0004 

TIFA0005 

TIFA0006 

TIFA0007 

TIFA00w8 

TIFA0009 

TIFAOOlO 

TIFAOOll 

T1FA0012 

TIFA0013 

TIFA0014 

TIFA0015 

TIFA0016 

TIFA0017 
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APPENDIX  II 
ILLUSTRATIVE  PROBLEM 
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Let's  consider  a  typical  mission  profile  for  a  tactical  aircraft.  This  mission 
will  cons  '  3t  of:  a  low  level  penetration  run,  a  pullup  over  the  target,  a  military 
power  climb,  a  roll  and  pullover,  a  dive  followed  by  a  5  "g"  pullup,  a  low  level 
dash  followed  by  a  4  "g"  pullup  into  a  military  power  climb  and  return  cruise 
to  a  normal  landing.  The  flight  description  is  taken  from  reference  9  and  is 
given  in  Table  II  and  loads  description  in  Table  III.  Assume  that  during  the  low 
level  penetration,  the  aircraft  sustains  damage  which  may  be  modeled  as  a 
through  crack  of  length  1. 74  inches.  The  input  cards  for  this  problem  are 
given  in  Table  IV  and  the  output  is  in  Table  V.  If  the  damage  is  modeled  as 
a  through  crack  of  length  2. 103,  the  input  changes  are  shown  in  Table  VI  and 
the  output  is  in  Table  VII.  Comparisons  with  analytical  results  from  reference 
9  show  that  both  analyses  predict  the  same  number  of  cycles  to  failure  for 
both  cases. 
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TABLE  II 

MISSION  PROFILE  DESCRIPTION 


Condition 


Description 

1 

Air  sp0t  d 

Low  Level  Dash 

450  kt 

180®  Turn 

450  kt 

Low  Level  Dosh 

450  kt 

180®  Turn 

450  kt 

Low  Level  Dash 

450  kt 

ISO®  Turn 

450  kt 

Low  Level  Dash 

450  kt 

180®  Turn 

450  kt 

Low  Level  Dosh 

450  kt 

180®  Turn 

450  kt 

Low  Level  Dosh 

450  kt 

180®  Turn 

450  kt 

Low  Level  Dosh 

4  50  kt 

Pullup 

450  kt 

Militory  Power  Climb 

0.87M 

Roll  ond  Pullover 

0.87  M 

Dive 

0.87  M 

Pullup 

0.87M 

Militory  Power  Climb 

0.87  M 

Cruise 

0.87M 

Descend 

0.87M 

Low  Level  Flight 

175  kt 

180®  Turn 

175  kt 

Low  Level  Flight 

175  kt 

180®  Turn 

175  kt 

Finol  Descent 

145  kt 

Numbir  Gust  Cycits 


An  2=0.26  An2  =  0.42 
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TABLE  m 


LOADS  FOR  MISSION  PROFILE 


Ao- 

P 

(p*i ) 

"  ( 

4  650 

0 

2  000 

0.568 

3  240 

0.  408 

6470 

0,256 

2  000 

0.  770 

3  240 

0.  652 

2000 

0.588 

3240 

0.  408 

6470 

0.  256 

2000 

0.  770 

3240 

0.  652 

2000 

0.  588 

3  240 

0.  408 

6  4  70 

0.256 

2000 

0.770 

3240 

0.  652 

2000 

0.  588 

3240 

0.4. Hi 

6470 

0.25ij 

2  000 

0.  770 

3240 

0.  652 

2  000 

0  588 

3  240 

0.408 

6  4  70 

0.256 

2  000 

0.770 

1  1 

2000 

0.5  88 

32  40 

0.408 

12 

6470 

0.  256 

12 

2000 

0.  770 

12 

3240 

0.652 

13 

2000 

0.  588 

13 

3240 

0.  408 

14 

131  70 

0.  145 

I5-I6 

1 1450 

0.  009 

17-18 

19150 

0. 00  5 

41 

2000 

0.  194 

41 

3240 

0.  031 

42 

5370 

0.  019 

42 

2000 

0.  5  88 

43 

2000 

0.  194 

43 

.3240 

0.  031 

44 

5370 

0.  019 

44 

2000 

0.  588 

44 

3240 

0.  408 

45 

7090 

0.240 

45 

2000 

0.  770 

45 

3240 

0 . 652 

46 

2000 

0.  588 

46 

in'40 

0.  408 

47 

7C90 

0 . 240 

47 

2000 

0. 770 

3240 

0.  652 

N 

(eyelta) 


I*' ) 


AFFDL-TR-70-107 


TABLE  IV 

DATA  DECK  FOR  PROBLEM  1 


»i4u't>t7U  mo  »«» 

I  I  i  M  I  I  I  I  I  I  I  i  I 


0TFDTRj‘UL65100 
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TABLE  IV  (CONTD) 


m 

a 

IS 

n 

HI 

a 

1 

3 

!•  t 

I 

ID 

1 

EQ 

IP 

jnur. 

iPi 

Em] 

■■ 

I] 

1 

n 

1 

21 

1 

II  at 

i, 

IT 

D 

[3 

n 

1 

m 

1 

IE 

1 

HE 

Rl 

IM 

D 

irai 

at 

40* 

m 

-1 

II 

II 

E 

l 

g 

II 

1 

K 

IIPB 

31 

1 

1 

1 

IP 

1 

1 

1 

|l 

1 

11 

II 

II 

H 

! 

1 

II 

1 

V 

\m 

si 

1 

1 

1 

K 

R 

1 

11 

1 

II 

3 

II 

II 

B 

1 

g 

II 

1 

IP 

IBE 

H 

1 

1 

1 

IP 

1 

1 

II 

II 

■ 

9 

II 

II 

1 

! 

i 

II 

1 

K 

m 

1 

1 

1 

1 

1! 

1 

1 

1 

n 

1 

1 

a 

II 

II 

N 

1 

1 

1 

1 

K 

IKE 

A 

1 

■P 

1 

1 

1 

II 

II 

H 

1 

II 

j 

g 

g 

1 

1 

i 

Wf 

31 

1 

1 

11 

i 

1 

R 

R 

III 

II 

II 

a 

a 

B 

1 

1 

1 

i 

1 

1 

m 

a 

1 

II 

1 

1 

1 

III 

PI 

1 

1 

P 

5 

g 

1 

1 

i 

Ini! 

Pi 

1 

1 

1) 

i 

I 

1 

1 

III 

a 

PI 

1 

1 

g 

1 

1 

¥ 

!IS* 

31 

1 

1 

r 

■1 

! 

1 

1 

n 

1 

1 

S 

1 

i 

1 

1 

1! 

9KE 

!B| 

1 

1 

V 

II 

1 

1 

1 

1 

1 

11 

II 

1 

PI 

a 

1 

E 

B 

1 

1 

i 

\m 

31 

1 

I 

i 

I 

1 

1 

I 

1 

II 

II 

1 

n 

P 

II 

1 

S 

g 

1 

1 

II 

im. 

1 

1 

» 

1 

1 

I 

1 

II 

1 

R 

II 

1 

1 

B 

II 

1 

i 

SKI 

SI 

1 

1 

IP 

I 

1 

1 

II 

1 

/ 

i 

11 

■ 

1 

I 

. 

II 

H 

9IB' 

n 

■k 

1. 

1 

1 

1 

!l 

1 

B 

PI 

i 

B 

3 

B 

I 

1 

11 

1! 

13! 

f'l 

/ 

R 

1 

! 

ill 

1 

5 

e 

B 

!l 

! 

I! 

IIP 

.si 

i 

1 

1 

I 

H 

5! 

M! 

PI 

R 

g 

II 

I 

i 

IB! 

± 

/ 

• 

! 

1 

1 

n 

1 

1 

1 

/ 

IS! 

[E 

I 

\ 

3 

11 

1! 

5  . 

>s 

1 

1 

I 

! 

II 

1 

1 

15 

1! 

11 

i 

IBS 

t . 

1. .. 

- 

A 

— 

- 

I 

!! 

! 

i 

E 

B 

I 

3 

B 

3 

II 

I! 

|i 

1 

i 

P 

P 

- 

1 

f 

ii 

I 

! 

I 

g 

11 

§ 

B 

1 

i> 

I| 

ligp 

.9 

I 

I 

1 

1 

ii 

q 

II 

B 

II 

i 

linr 

5! 

1 

1 

i 

• 

1 

11 

!i 

5 

s 

II 

1 

1! 

IIP 

mt 

1 

1 

.J 

! 

1 

1 

1 

5 

n 

e 

i 

U 

I 

1 

Ibp 

1/ 

1 

1 

ii 

1 

I 

11 

i 
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TABLE  IV  (CONTD) 


o  it  a  s  4 'I  |<  IT  (  •  10  ti  ti  tauta'utait  laae  tiaaaaa4aa  ataraa  aaae  a<aa  saa4aaa4a7a(a*  40. 

o  - - 


TABLE 


AFFDL-TR-70-107 


o 

I 

N 

(/) 

4/1 


CC 

i/t 

Q 


> 


VA 

o 

Ul 

O' 

O' 

o 

O' 

O' 

O' 

«■» 

0 

0 

• 

< 

o 

H 

UJ 

0 

n 

• 

• 

u 

• 

0 

# 

o 

A 

4/1 

N 

X 

o 

< 

z 

4b 

4b 

D 

0* 

0 

Z 

O' 

* 

M  O 

O' 

z  o 

o 

3  • 

O' 

«J  0 

0 

< 

N 

< 

M 

K 

UJ 

0 

• 

0 

z 

H  Z 

o 

4b 

« 

o 

1  «J 

o 

in  .j 

< 

V) 

N  < 

H 

z 

# 

4/> 

o  z 

o 

*>« 

M 

N  in 

z 

u. 

a 

X 

O 

o 

oc 

Z 

0 

z 

UJ 

UJ 

in 

a 

0 

3 

in 

z 

<N 

# 

o 

4/1  ^ 

z 

in 

0 

H  1 

4J 

ac 

h* 

• 

Z  4b 

< 

a 

CC 

<  O 

z 

h" 

< 

UJ 

H  O 

o 

o 

4/)  O 

< 

X 

‘  z  o 

Ur 

u. 

lb 

o 

o  o 

O 

Wirt 

O  0 

< 

Z 

u. 

z 

o 

N 

^  '/I 

-J 

<  • 

-J 

Z 

< 

M  O 

< 

u 

o 

QC 

M 

UJ 

< 

M 

UJ 

H 

oc 

& 

W  tl 

oc 

>- 

< 

Z 

a 

4b 

H 

z  u 

u 

o 

I 

rg 

< 

►- 

Ui 

CC 
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FINITE  WIDTH  CORRECTION  FACTOR 

BETA(2)  *  S0RT(2/(PI*A/8}*TAN(PI*(A/3)/2)| 

WHERE. ...B  >  0.80000000E  01 


IW  ‘ 
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